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Abstract 

We describe the method for extracting the elastic scattering phase shift from a lattice 
simulation at an introductory level, for non-lattice practitioners. We consider the scattering 
in a resonant channel, where the resulting phase shift 5(s) allows the lattice determination 
of the mass and the width of the resonance from a Breit-Wigner type fit. We present the 
method for the example of P-wave 7T7T scattering in the p meson channel. 

1 Introduction 

The determination of the strong decay width of a hadronic resonance in lattice QCD is 
a much more demanding task than the determination of its approximate mass. The only 
available method (that was applied up to now) was proposed by Liischer [1] and is rather 
indirect. It applies for the case when the resonance appears in the elastic scattering of two 
hadrons HiH 2 — > R — > H 1 H 2 . 

• First, the energy spectrum E n of the system of two interacting hadrons H\H 2 enclosed 
in a few-fermi box has to be determined. The system is illustrated in Fig. [TJ The 
spectrum in a finite box E n is discrete and few (one or two) lowest energy levels have 
to be determined by lattice simulation. 

• The shift of the energy E n with respect to the non- interacting energy i?m(Pi) + 
Eh2(P2) (E Hi (pi) = \J mf + pi 2 ) gives info on the interaction between Hi and H 2 . 
Liischer derived a rigorous relation between the energy shift E n — Em — Eh2 and the 
elastic phase shift 8(s) for H\H 2 scattering in continuum pp. The measured energies 
E n can be used to extract the phase shift 5(s) evaluated at s = E% — P 2 , where E n is 

1 Presented by S. Prelovsek at Bled Mini- Workshop 2011, Understanding Hadronic spectra, 3-10 July, 
Bled, Slovenia. 
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Figure 1: The energy of two hadrons in a box of size L. On the left, L > fm and E(L) ~ 
E m (p 1 )+E H2 (p2)- On the right, L ~ few fm and energy gets shifted due to their interaction, 
i.e. E(L) ~ E H1 ( Pl ) + ^ 2 (p 2 ) + AE(L). 



the energy of the system and P its total momentum. In order to extract S(s) at several 
different values of s, the simulations are done for several choices of total momenta P 
of the H\Hi system, which leads to different values of s = E\ — P 2 . 

• The resulting dependence of S(s) as a function of s can be used to extract the mass rriR 
and the width Tr of the resonance R, which appears in the elastic channel H1H2 — » 
R — > HiH 2 . For this purpose, the S(s) can be fitted with a Breit-Wigner form or some 
other phenomenologically inspired form, which depend on wir and r^. 

The described method, needed for the determination of the resonance width Tr, is rather 
challenging. It requires very accurate determination of a few lowest energy levels of the 
system HiH 2 , since the resulting phase shift depends ultimately on the energy shift. Among 
all the meson resonances, this method has been up to now rigorously applied only to p 
resonance. Although Liischer proposed the method already in late 80's [Tj, the first lattice 
attempt to employ it to hadronic resonances had to wait until 2007 J2]. Since then, several 
studies of p have been carried out [31 H], with the most up to date ones [5l [7J. 

This talk briefly describes the method to extract S(s), and Tr on an example of tttt 
scattering in the p channel. It is based on a recent simulation [6], which is the statistically 
most accurate determination of any strong meson width on one lattice ensemble. The purpose 
of this talk is to highlight the main physical reasoning, which lies behind the lattice extraction 
of S(s), rriR and Tr, omitting most of technical details. 

The sections follow the order of steps required, which are listed as items in the in- 
troduction. Section II describes the determination of spectrum E n of the coupled system 
HiH 2 -H- R. The Section III described why E n allow one to extract the elastic phase shift 
S(s). The extraction of the resonance parameters rriR and Tr from the phase shift 5{s) is 
done in Section IV. We end with conclusions. 
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2 Spectrum of two hadrons in a finite box 



The p meson is a resonance in 7T7t scattering in P-wave, and has quantum numbers I G (J PC ) = 
1 + (1 ). The total momentum P of the coupled ttit — p system can have values -j^d , d G Z 3 
due to the periodic boundary condition in the spatial direction, and we use the following 
three choices 

P = (0, 0, 0) , ^-(0, 0, 1) , fj-(l, 1, 0) and permutations . (1) 

This enables us to obtain several values of s = — P 2 for the system, thereby allowing the 
determination of 5(s) for these values of s without changing the spatial volume. 

Our simulation is performed on an ensemble of 280 [8J gauge configurations with dy- 
namical u/d quarks, where the valence and dynamical quarks employ improved Wilson- 
Clover action. The corresponding pion mass is m n a = 0.1673 ± 0.0016 or m v = 266 ± 4 
MeV. The lattice spacing is a = 0.1239 ± 0.0013 fm and we employ a rather small vol- 
ume iVf x Nt = 16 3 x 32, which allows us to use the costly full distillation method [9] for 
evaluating the quark contractions. 

On the lattice, the discrete energies of the system E n can be extracted after computing 
the dependence of the correlation matrix Cij(tf,ti) on Euclidean time tf — U 

C^U) = (O|0 4 (f,) 0){UM = £<WH<9}> e-^'/-*) . (2) 

n 

The analytical expression on the right is obtained by inserting the complete set Yl n \ n )( n \ 
of physical states n with given quantum numbers. The interpolators Oi have the quantum 
numbers of the system in question. In our case the interpolators have quantum numbers 
J PC = 1 and | J, -Z3) = 1 1,0) and total three-momentum P. They have to couple well to 
the 7T7T state and the quark-antiquark resonance p. 

For each choice of P (pQ), we use 16 interpolators, listed in detail in Eq. (21) of [6]. We 
employ fifteen interpolators of quark-antiquark type 

Of{t) = ^ 7f[^ M (^ x ) + d ^d (t,x)] , (3) 

X 

where Ti denotes different color-spin-space structures with the same resulting quantum num- 
ber J PC = 1 and I J, I3) = |1,0). We use also one vr(pi)7r(p 2 ) interpolator, where each 
pion is projected to a definite momentum 

° n7T (t) = ^k + (Pi)vr-(p 2 )-7r-( Pl )7r + (p 2 )] , p 1+ p 2 = P ^(pi) = J>^ x 975^ (t, x) 

X 

(4) 

In practice, the interpolator is the most important among our 16 interpolators, since it 
couples to the scattering state much better than the quark-antiquark interpolators. Let us 
note that all other lattice studies aimed at T p used at most one quark-antiquark and one irn 
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Figure 2: Contractions for 1=1 correlators with qq (j3]) and ittt (jj]) interpolators. 

interpolator, which may not always allow for reliable extraction of the first excited energy 
level E 2 . 

Given the 16 interpolators, we compute the 16 x 16 correlation matrix Cij(tf,U) for 
all initial and final time-slices = 1,..,Nt = 32. The needed Wick contractions that 
enter the correlation matrix with our qq and nir interpolators are depicted in Fig. |2j The 
contributions (a,c,e) in Fig. [2] cannot be evaluated solely from the quark propagator from 
one point (tj,Xj) to all other points of the lattice (such a propagator allowed most of the 
spectroscopy studies in the past). The contributions (a,c,e) require the propagators from all 
and to all points on the lattice, which is too costly to evaluate in practice. We use the recently 
proposed distillation method for this purpose j9], which enables the exact computation of 
the required contractions. 

We average the resulting correlators (i) over all initial time slices £j at fixed time sepa- 
ration tf — ti, (ii) over all directions of momenta P ([Q) and (iii) over all directions of the p 
meson polarization. 

The time dependence tf — ti of the correlators Cij(tf,ti) (|2j) contains the information on 
the energies of the system E n , and several methods for extracting E n from are available. 
We extract two lowest energy levels E n=12 of the system from the 16 x 16 correlation matrix 
Cij(tf,ti) using the so called variational method [10J, which is the most established among 
the available methods. Table [1] displays the extracted lowest two energies E n= i t2 of the 
coupled 7T7T — p system for our three choices of total momenta P flTJ . 

The spectrum E n in Table CD for our finite box is the main result of this section. Each 
energy level corresponds to a different value of s = E\ — P 2 , as calculated from E n and P 
in the Table [TJ In fact, the table lists values of s obtained from the discrete lattice version 
of the dispersion relation, which takes into account part of the corrections to s = E\ — P 2 
due to finite lattice spacing [6]. 

3 Extraction of the phase shifts from energy levels 

Let us consider the case when the resonance R can strongly decay only to two spinless 
hadrons H\ and H 2 , so one has elastic scattering of H\ and H 2 . We point out that the 
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E n a 


sa 
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f (0,0,0) 


1 


0.5107(40) 


0.2608(41) 


130.56(1.37) 


f (0,0,0) 


2 


0.9002(101) 


0.8103(182) 


146.03 (6.58) [*] 


■f (0,0,1) 


1 


0.5517(26) 


0.1579(29) 


3.06 (0.06) 


f (0,0,1) 


2 


0.6845(49) 


0.3260(69) 


156.41(1.56) 


f (1,1,0) 


1 


0.6933(33) 


0.1926(49) 


6.87(0.38) 


1(1,1,0) 


2 


0.7868(116) 


0.3375(191) 


164.25(3.53) 



Table 1: The results for two lowest levels n — 1, 2 of the coupled tttt — p system with three 
choices of total momentum P on our lattice with m n a = 0.1673 ± 0.0016, L = 16a and the 
lattice spacing a = 0.1239 ± 0.0013 fm. The energy levels E n are obtained by multiplying 
E n a with a -1 ~ 1.6 GeV. The invariant mass squared of the system is s = E 2 — P 2 , but the 
dimensionless value in the table s a 2 is obtained using the discretized version of this relation 



non-elastic case, when a resonance can decay strongly to several final states (i.e. HiH 2 and 
H[H' 2 ), is much more challenging for a lattice study. 

Suppose one encloses two hadrons £fx(pi) H 2 (p2) with three-momenta pi and p 2 into a 
large box of size L > fm and measures their energy. In a large box, they hardly interact 
and their energy is equal to sum of individual energies E non ~ mt = Ehi{px) + Ehi^P?) with 
Eh(p) = \J m2 H + P 2 - Now, let's force -ffi(pi) and H 2 (p2) to interact by decreasing the size of 
the box to L of a few fm. The energy of the system E(L) = Euiipi) + Eh2{j>2) + AE(L) is 
shifted with respect to E non ~ mt : it will increase (AE(L) > 0) if the interaction is repulsive 
and decrease (AE(L) < 0) if the interaction is attractive. This simple physical reasoning 
indicates that the energy shift AE(L) gives info on the interaction. 




Figure 3: The scattering of two interacting particles as series of the interaction vertex M(Sl) 
and the scattering of non- interacting particles F at finite L [11]. 



In fact, the energy shift AE(L) and the energy itself E(L) do not only give us "some" 
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info on the interaction. According to the seminal analytic work of Luscher [1], E(L) or 
AE(L) rigorously tells us the value of the elastic scattering phase shift of HiH 2 scattering 
at L — > oo, i.e. 8(L = oo): 

Luscher method : E(L) — > S{s, L = oo) s = E(L) 2 - P 2 (5) 

The derivation and the resulting formulae between E(L) and 5 are lengthy and rather com- 
plicated, but let us briefly explain at least why E(L) contains info on S(L = oo). A nice and 
clear quantum-filed theory derivation is given in [TT] and the main message is illustrated in 
Fig. El The scattering of two interacting spin-less hadrons HiH 2 at finite L (for degenerate 
case rriHi = tuh2 = Tn) is represented in QFT by series of: 

• scattering of two non-interacting hadrons at finite L, represented by F. The expression 
F contains sums over the loop momenta k, which are allowed in a finite box L with 
periodic boundary conditions in space. Here /(A; ,k) stands for dependence of the 
vertices on the left and right on ko and k. 

• the interaction vertex M with four hadron legs. This vertex depends on the elastic 
phase shift Si (at infinite volume) for the case of elastic scattering in the Z-th partial 
wave. 

The physical scattering requires resummation of the bubbles in Fig. [3j with non-interacting 
parts F and the interacting parts M, giving AF —jrp A'. The positions of the poles of the 
sum AFj^jpA' obviously depend on M and therefore on Si. The positions of the poles 
dictate the possible energy levels of the system E n (L), so the energy levels E n (L) depend on 
M and therefore on Si. 

The purpose of the above illustration was just to indicate why E n (L) depend on Si. In 
the case of tttt with J p = 1~, the relevant wave has I = 1 and we denote the corresponding 
phase by S = Si. The complete analytic relations between E n (L) and S(s) needed for our 
case of the tttt scattering with J PC = 1 and 1 = 1 are provided in [6] (for every |P| a 
different form of relation applies). These allow to extract S for each of our six energy levels 
in Table [Hand the resulting phase shifts are given in the same Table. 

The presented Luscher formalism applies only for the case of elastic scattering. The tttt 
state is the only scattering state in this channel for energies when 4tt state cannot be created, 
i.e., when s = E 2 < (4m 7r ) 2 . For our m^a = 0.1673 this is valid for all six levels, with the 
exception of the level E\ at P = 0, which is above 4tt inelastic threshold. As the Luscher 
analysis is not valid above the inelastic threshold, we omit this level from further analysis. 

The resulting scattering phase shifts for five values of s are shown in Fig. HI This is the 
main result of the lattice study; the resonance properties will be obtained by fitting S(s) in 
the next section. 

Note that the resulting phases are determined with a relatively good precision, which 
is better than in other available lattice studies of p at comparable u/d quark masses. The 
good precision can be traced back to various advanced techniques we used: the distillation 
method for evaluating contractions, usage of a large interpolator basis and average over all 
initial time slices, directions of momenta P and polarizations of p. 
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Figure 4: The nir phase shift S(s) (in degrees) for five different values of dimensionless 
so 2 = (E n a) 2 — (Pa) 2 , extracted from our lattice study [6J. The s is obtained by multiplying 
so 2 with (a" 1 ) 2 ~ (1.6 GeVf. 



4 Extracting resonance mass and width from the phase 
shift 

The phase shift 5{s) in Fig. HI obtained directly from the lattice study, can be used to extract 
the properties of the resonance, in our case the p. The phase shift has a typical resonance 
shape: it passes from 5 ~ 0° to 5 ~ 180°: the point where it crosses 90° gives the position of 
the resonance (s = m 2 ), while the steepness of the rise gives its width T p . In particular, 5 is 
related to resonance parameters by expressing the scattering amplitude a\ in terms of 5 on 
one hand, and with Breit-Wigner form in the vicinity of the resonance on the other hand 

-v^r(s) e 2lS ^ - 1 , < 

ai = T^JTv^nY) = 2% ' (6) 

Relation ((H]) can be conveniently re-written as 

y/sF(s) cot 6(s) = m 2 — s . (7) 

The decay width significantly depends on the phase space and therefore on m n , so the decay 
width extracted at m n ~ 266 MeV could not be directly compared to the measured width. 
So, it is customary to extract the p — > iyty coupling g p7T7T instead of the width, where the 
width 



P * 3 ^ T n = T(m 2 " 



s bn 

depends on the phase space for a P-wave decay and the coupling g p7T7T . The coupling is 
expected to be only mildly dependent on m n , which was explicitly confirmed in the lattice 
studies [5j [7] and analytic study [12] . In flSJ), p* denotes the pion momentum in the center- 
of-momentum frame and we extract it from s using a discretized version of relation yfs = 
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Figure 5: The crosses are the inr phase shift 5(s) (in degrees) for five different values of 
dimensionless sa 2 = (E n a) 2 — (Pa) 2 , extracted from our lattice study [6]. The line is the 
Breit- Wigner fit fl7|8|) for the resulting m p and g pn7T in Table [2j The physical value of s is 
obtained by multiplying sa 2 with (a -1 ) 2 — (1.6 GeV) 2 . 





lattice (this work [6]) 


exp [PDG] 




m w ~ 266 MeV 






792 ± 12 MeV 


775 MeV 




5.13 ±0.20 


5.97 



Table 2: Our lattice results for the resonance parameters [5], compared to the experimental 
values. 

2^/m 2 + p* 2 [6j. Inserting T(s) §E§ into ([7|), one obtains an expression for 5(s) in terms of 
two unknown parameters: m p and g pn7T . We fit these two parameters using five values of S(s) 
given in Fig. H] and Table [IJ and we get the values of resonance parameters in Table [2] with 
small statistical errors. 

The resulting p-meson mass in Table [2] is slightly higher than in experiment, as expected 
due to TYi-ff = 266 MeV > nf* v . The coupling g p7T7T is rather close to the value gf%* derived 
from the experimental width Y p xp . 

5 Comparison to other lattice and analytical studies 

The comparison of our results for m p and Y p to two recent lattice studies [7] is compiled 
in Fig. 8 of [7]. Our result has the smallest error on a given ensemble, demonstrating that 
accurate lattice determination niR and Tr for (some) resonances is possible now. The other 
two lattice studies are done for two [7] and four [5] pion masses and explicitly demonstrate 
mild dependence of g p7T7T on m w . The discussion concerning the (dis) agreement of the three 
lattice studies is given in [7] and will be extended in [T5] . 
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The comparison of our S(s) to the prediction of the lowest non-trivial order of unitarized 
Chiral Perturbation Theory [T3] is given by the solid line in El which has been recalculated 
for our m n = 266 MeV in [T5] . The lowesto order prediction does not depend on unknown 
LECs and agrees reasonably well with our lattice result, given by the bullets. 
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Figure 6: The irn phase shift in the p channel Su(p) = 5{p*) at m n = 266 MeV: the solid 
line (indicated by "Unit 0(p 4 )") gives prediction of the lowest order of Unitarized Chiral 
Perturbation Theory [HI [15] , while bullets are our lattice data. 



6 Conclusions 

We highlighted the main physical reasoning, which lies behind the lattice extraction of elastic 
phase shifts 8(s) and the resonance parameters rriR and Y R . The purpose was to present the 
general principle of the method and omit the technical details. The method was presented 
on the example of titt — > p — > tttt scattering. This example demonstrates that a proper 
first-principle treatment of some hadronic resonances on the lattice is now possible. 
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